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Abstract 

Using  an  automatic  relay  computer,  a  method  for  determining  the  inter¬ 
ference  color  of  thin  films  is  studied.  Calculation  of  chromaticity  of  the 
interference  color  is  performed  more  conveniently  by  the  summation  of  inte¬ 
grand  at  equal  intervals  of  wavenumber  than  at  equal  intervals  of  wavelength, 
if  the  change  of  refractive  index  with  wavelength  is  neglected. 

The  CIE  colorimetric  data  of  Efxx,  Exi/x,  ExZx  are  transformed  into  a 
table  using  the  wavenumber  in  200  cm~*  intervals  as  variable  and  are  given 
in  the  range  extending  from  13,000  cm~^  to  27,000  cm'L 

1.  Introduction 

The  colorimetric  calculation,  needed  to  determine  the  chromaticity  of  a  color 
from  its  spectral  intensity  distribution,  is  rather  laborious.  Mechanical  computers 
suited  to  this  purpose  have  been  designed  and  used.  Computation  of  the  inter¬ 
ference  color  of  thin  films  is  much  more  tedious  than  that  of  ordinary  color.  A 
method  to  save  the  labor  and  simplify  the  theoretical  interpretation  of  results  has 
been  adopted  by  neglecting  the  effect  of  multiply  reflected  rays“*'. 

Lately,  such  high-speed  automatic  calculating  machines  as  the  electronic 
computer  and  the  electric  relay  computer  have  been  developed.  It  seems  natural 
that  much  time  and  labor  of  above-mentioned  computation  will  be  saved  by  using' 
these  computers.  An  application  of  IBM  electronic  computer  to  colorimetric 
•calculation  has  been  reported*’. 

Pressed  by  the  need  of  determining  the  interference  color  of  stained  glass 
surfaces  coated  with  anti-reflecting  Aims,  an  attempt  has  been  made  to  apply  an 
•electric  relay  computer  to  this  problem. 

The  computer  FACOM-128,  made  by  the  Fuji  Communication  Apparatus 
MFC.  Co.,  was  used  simply  because  it  was  easily  accessible  to  the  author,  and 

1)  H.  KuboU:  J.  Opt.  Soc.  Am.  40  (1950)  146. 

2)  H.  Kubota:  Rep.  Inst.  Indust.  Sci.  Tokyo  Univ.,  2  (1952)  205. 

G.  C.  Monch:  Optik  9  (1952)  97. 

3)  W.  E.  White  and  D.  L.  MacAdam:  J.  Opt.  Soc.  Am.  47  (1957)  605. 
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after  all,  it  was  found  for  the  present  most  convenient  as  it  has  a  large  capacity 
of  memorial  elements.  For  calculation,  the  program  of  colorimetric  procedure 
was  devised  to  be  easily  applicable  to  the  operation  of  the  relay  computer. 

2.  The  Colorimetric  Calculation  Method 

Description  of  the  procedure  adopted  by  the  CIE  Committee  on  Colorimetry 
to  specify  the  chromaticity  of  a  color^’  is  here  left  out  except  on  a  few  points 
considered  necessary  for  explanation. 

Let  Ek  be  the  energy  per  unit  width  of  wavelength  in  the  distribution  of  the 
light  source  energy,  Rx  the  reflection  coefficient  of  the  surface  under  consideration 
for  the  light  of  wavelength  X,  Xx,  Vx,  zx  the  tristimulus  values  of  the  spectrum 
color,  X,  Y,  Z  the  tristimulus  values  of  the  reflected  light,  x,  y,  z  the  tristimulus 
coefficients  of  the  reflected  light,  then  the  chromaticity  of  the  reflected  light  is 
given  as 

A'= ,  Y=^ExVxRxdX,  Z=^ExZxRxdX  (1) 

V=XIS,  y=YIS,  z=Z/S 

S=X+Y+Z  (2) 

Usually,  numerical  calculations  by  summation  are  substituted  for  the  integra¬ 
tions  in  (1).  Although  there  are  two  ways  to  perform  this  calculation,  that  is,  the 
selected  ordinate  method  and  the  weighted  ordinate  method,  the  latter  was  used 
by  the  reason  described  below. 

The  reflection  coefficient  Rx  of  a  double  layer  is  given  by 

=  1  -  (1  -  r,*)(l  -  ri*)(l  -  r,*)/Z) 

+ 2  { r#*)  cos  dz  +  »'on(l  +  >'**)  cos 

-I-  r^rdcos  (di  -1-5*) + n*  cos  (5i  -5*)]}  ( 3 ) 

where  r,,  r,,  r*  are  the  amplitude  reflection  coefficients  of  the  boundary  planes  between 
mediums  as  shown  in  Fig.  L 

For  the  sake  of  simplicity,  only  the  case  of  normal  incidence  will  be  considered 
hereafter.  Let  the  refractive  indices  of  the  mediums  be  tto,  tti,  n*,  n,  in  ascending 
order  from  the  bottom,  we  have  the  following  relations, 

4)  Committee  on  Colorimetry;  J.  Opt.  Soc.  Am.  34  (1944)  633. 

A.  C.  Hardy:  Handbook  of  Colorimetry.  (The  Technology  Press,  M.I.T.,  Cambridge, 
1936). 
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Class  n. 

Fig.  1.  Structure  of  a  double  layer 


•  «o+«i  ’  *  «!  +  «*’  *  «2  +  «S 


(4) 


The  phase  angles  in  (3)  are  given  as 

d,  =  {AnlX)n,di,  9^={4n|X)n■^d^,  *  (5) 

where  di,  dt  are  the  thicknesses  of  the  mediums  1  and  2.  Let  us  use  m/x  as  the 
imit  of  measure  for  wavelength  and  optical  thickness  of  Alms. 

If  the  interference  color  of  this  doubly  layered  surface  is  to  be  determined, 
firstly  the  reflection  coefficient  Rx  has  to  be  calculated  for  the  wavelengths,  which 
are  selected  as  the  calculation  points.  Then,  in  the  case  of  selected  ordinate 
method,  Rxs  must  be  summed  up  for  all  the  selected  wavelengths  to  obtain  the 
tristimulus  values  X,  Y,  Z.  In  the  case  of  weighted  ordinate  method,  the  sum¬ 
mation  is  to  be  performed  after  multiplying  Rxs  with  ExXx,  ExVx,  ExZx  of  cor¬ 
responding  wavelengths  as  weighting  factors.  Rx  involves  £i,  3t  in  the  forms  of 
cos^i,  cosdf  C06(di-l-<x).  cos(di—8t)  as  the  expressions  (3)  show.  When  the 
computation  is  made  with  a  hand  operated  calculating  machine,  numerical  values, 
of  these  functions  are  found  in  a  trigonometric  function  table.  But  for  the 
automatic  calculating  machine,  they  will  have  to  be  calculated  separately.  Using 
a  calculator  such  as  a  relay  computer,  the  speed  of  which  is  rather  slow,  not  only 
it  makes  the  programming  troublesome  but  also  the  time  of  computation  becomes 
so  long  that  the  advantage  of  using  such  an  automatic  computer  will  be  counter¬ 
vailed. 

In  the  selected  ordinate  method,  the  wavelengths  of  the  calculation  points  are 
given  in  the  intervals  which  are  inversely  proportional  to  the  values  of  ExXx, 
ExVx,  ExZx  and  therefore  there  is  no  simple  rule  to  define  their  intervals.  In  the 
weighted  ordinate  method,  they  are  given  in  equal  intervals,  for  example,  at  every 
lOm/i.  As  it  is  apparent  from  (5),  8i  and  8^  do  not  change  in  equal  intervals 
even  if  A  changes  with  equal  intervals.  Neglecting  the  change  of  refractive  indices 
iti,  with  the  wavelength.  Si  and  8t  become  proportional  to  1/i. 

Therefore,  if  the  quantity,  which  changes  proportionally  to  1/>1,  for  example. 
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the  wavenumber  v=CiA  (where  C  is  equal  to  10“  m/i),  is  taken  as  a  variable  and 
changes  with  equal  intervals  of  Jv,  both  S,  and  change  with  equal  intervals  of 
Aoi  and  respectively,  where 

Jdi  =  {4n/C)Ji>'nidi,  (6) 

Hence,  if  sin  5,,  sin  Jd,,  cos  5,  and  cosi<5i  are  given,  sin(5i  +  J«)i)  and  cos(5,+J5i) 
are  determined  by  using  the  sum-rule  of  trigonometrical  founctions 

sin  (51+ J5i)=sin  5i  cosi5i+cos5,  sin  45,  x 

L  (7) 

cos  (5i+ J5,)=cos5i  cos  J5,— sm  5i  sm  J5i  ’ 

Repeating  this  process,  sin(5i+mJ5,)  and  cos(5i+mJ5i)  are  determined  easily, 
where  m  is  an  integer. 

As  this  process  is  of  only  arithmetical  operations,  there  is  no  trouble  in  using 
an  automatic  computer,  and  the  computation  was  performed  on  an  automatic 
relay  computer  by  the  weighted  ordinate  method  taking  the  wavenumber  in  equal 
intervals  as  the  variable. 


3.  The  Computation 

The  CIE  colorimetric  data  give  the  numerical  values  of  ExXk,  E^y^,  ExZx  for 
the  CIE  A-,  B-  and  C-Illuminants  in  lOm/i  intervals.  It  is  necessary  to  convert 
them  into  a  form  of  tables  with  equal  intervals  of  v.  At  first,  the  values  of  ExXx 
etc.  corresponding  to  the  wavelength  i=ClvmM,  where  v  is  a  tabulation  point, 
were  determined  applying  the  interpolation  method. 

As  the  interpolation  method,  the  third-difference  osculatory  method  was  used. 
This  method  was  first  used  by  Judd“’  to  interpolate  the  visibility  function  and 
thereafter  by  others*’  to  interpolate  the  values  of  ExXx  etc.,  given  for  every  10  m^, 
in  1  mu  intervals.  In  this  metivod,  the  function  in  the  interval  between  two 
successive  values  given  in  the  table  is  expressed  by  such  a  third  order  curve  that 
the  tangent  at  one  end  of  the  curve  is  also  the  tangent  to  the  curve  in  the  adjoining 
interval. 

Assume  that  the  values  of  a  function  of  >1,  V\,  are  given  in  equal  intervals  of 
lOm/i,  and  the  first,  second  and  third-order  differences  are  as  given  in  Table  1. 
Then,  Vx  for  -lo<^<>lo+10  is  given  by 

Vx=  Va,_io  +  A’iJi  Vx.-io+A'zJjVaj-io  +  A'jJ,  Va,_,o  (8 ) 

5)  D.  B.  Judd:  Bur.  Stand.  J.  Res.  6  (1931)  465. 

6)  A.  C.  Hardy:  loc.  cit.  p.  32. 
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where 


/Ci=(>l-i,-t-10)/10 

1 

10)(>l  -  4)/200 

l 

(• 

(9) 

TaWe  1. 

Vx  d^Vx 

diVx 

d»Vx 

10 

VXg_10 

JiVaq-io 

V*. 

dxVx, 

Vxg+20 

^2V^X. 

daVx^-iQ 

Using  this  equation,  the  values  of  EkXa.  etc.  for  the  desired  v  are  determind, 
but  these  are  the  values  per  unit  width  of  the  wavelength.  The  values  E-^x^  etc. 
needed  in  the  integration  by  v,  are  obtained  by  converting  them  to  the  values  per 
unit  width  of  the  wavenumber  by  multiplying  dXjdv.  That  is, 

E^^=ExXx^J^  =  —^,ExXx.  (10) 

dv  ir 

As  mentioned  above,  the  CIE  data  are  given  in  every  10  m/^.  It  does  not  tell 
where  the  value  becomes  zero  within  the  interval  of  lOm/i.  Therefore,  it  was 
assumed  that  the  value  -becomes  zero  at  a  wavelength  10  m^  away  from  the 
wavelength  at  which  the  value  is  the  nearest  to  zero.  In  the  CIE  data,  the 
wavelength  range  in  which  any  of  the  values  Efxx,  ExVx,  Efzx  does  not  become 
zero,  is  from  380  m/«  to  770  m/i.  Accordingly,  by  extending  lOm/i  at  both  ends, 
one  may  assume  that  the  values  are  given  in  the  range  from  370  m//  to  780  m/x. 
The  wavenumbers  corresponding  to  these  limits  are  27, 027cm"‘  and  12, 820  cm"*. 
The  range  of  wavenumber  in  which  the  values  of  £vXv  etc.  should  be  tabulated,* 
was  taken  as  from  13,000  cm"*  to  27,000  cm"*.  The  interval  in  wavenumber  was 
selected  as  200  cm"*,  corresponding  to  3m/i  in  the  short  wavelength  end  and 
to  12m/i  in  the  long  wavelength  end  so  that  the  corresponding  interval  in  wave¬ 
length  should  not  exced  lOm/ti  except  at  the  long  wavelength  end,  where  it  is 
usually  not  of  importance. 

The  conversion  was  performed  using  the  relay  computer.  The  results  for  the 
CIE  A-  and  C-Illuminants  are  tabulated  in  Table  2.  The  values  in  the  table  are 
multiplied  by  a  coefficient  Jv=200  for  convenience.  If  the  values  in  the  table  are 
summed  up,  they  give  the  tristimulus  values  X,  Y,  Z  of  the  illuminant.  If  the 
value  per  unit  width  of  wavenumber  is  needed,  they  must  be  divided  by  200. 

Directions  for  the  use  of  tables  are  as  follows.  Let  the  r  coefficient  eflection 
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Iluminant  A 

Illuminant  C  | 

cin”' 

Eax 

EaV 

Eaz 

Ecx 

Ecy 

Eci 

13000 

2.5 

0.1 

1.2 

0.1 

13200 

5.0 

2.4 

1.3 

1.2 

1 

13400 

9.2 

3.6 

2.9 

1.4 

i 

13600 

20.9 

8.2 

5.8 

2.5 

1 

13800 

43.1 

14.4 

13.5 

4.4 

1 

14000 

84.5 

30.0 

28.9 

10.0 

14200 

156-7 

56.5 

58.9 

20.9 

14400 

TXJA 

103.5 

117.1 

42.1 

14600 

542.9 

196.6 

239.9 

87.0 

14800 

946.9 

344.9 

446.6 

162.5 

15000 

1584.1 

581.3 

791.1 

290.4 

15200 

2548.4 

946.0 

1333.8 

495.5 

15400 

3765.6 

1424.6 

2045.6 

773.2 

15600 

5193.5 

2019.2 

2911.1 

1131.8 

15800 

6626.3 

2682.7 

3861.4 

1563.1 

16000 

8005.4 

3420.8 

1.1 

4853.2 

2073.6 

0.8 

16200 

9021.4 

4136.9 

2.5 

5688.4 

2608.9 

1.5 

16400 

9432.5 

4748.8 

3.0 

6207.4 

3124.9 

1.5 

16600 

9293.0 

5280.8 

6.3 

6425.8 

1  3648.5 

4.2 

16800 

8707.3 

5705.0 

7.8 

6405.8 

4199.2 

5.9 

17000 

7819.0 

6010.8 

8.9 

6177.7 

;  4752.5 

7.5  i 

17200 

6760.1 

6191.4 

11.7 

5765.6 

1  5281.2 

10.3 

17400 

5677.4 

6213.9 

12.5 

5213.3 

5711.1 

11.4 

17600 

4640.5 

6100.3 

14.6 

4573.2 

6011.6 

14.3 

17800 

3700.4 

5867.8 

20.0 

3884.1 

6163.9 

21.3 

18000 

2889.3 

5541.5 

30.3 

3194.4 

:  6124.8 

33.7 

18200 

2200.7 

5146.1 

47.4 

2534.1 

!  5924.3 

54.6 

18400 

1635.6 

4709.4 

74.1 

1940.2 

5583.9 

87.7 

18600 

1170.3 

4234.5 

110.5 

1424.2 

5147.3 

134.2 

18800 

792.4 

3728.2 

156.0 

985.4 

4634.8 

193.3 

19000 

488.9 

3200.4 

207.1 

626.7 

1  4109.3 

264.3 

19200 

257.7 

2661.7 

270.6 

346.0 

3575.8 

362.2 

19400 

111.9 

2120.7 

374.9 

158.7 

3041.0 

539.7 

19600 

31.8 

1619.4 

499.2 

49.2 

2532.5 

780.0 

19800 

5.5 

1216.2 

618.9 

10.6 

2085.4 

1063.9 

20000 

13.5 

896'.  0 

754.5 

26.0 

1700.5 

1432.0 

20200 

36.4 

668.9 

913.8 

76.9 

1392.7 

1910.7 

20400 

75.0 

504.9 

1106.9 

169.5 

1141.6 

2503.7 

20600 

126.6 

385.1 

1357.2 

307.3 

937.7 

3294.3 

20800 

186.9 

299.6 

1631.0 

480.6 

771.2 

4202.4 

21000 

253.7 

229.8 

1914.3 

692.5 

1  627.1 

5223.1 

21200 

321.3 

177.1 

2183.3 

923.0 

1  507.1 

6269.9 

21400 

373.7 

136.0 

2358.1 

1129.7 

,  411.5 

7122.9 

21600 

414.6 

106.0 

2460.9 

1320.6 

337.6 

7833.0 

21800 

434.6 

82.0 

2464.6 

1457.3 

1  276.0 

8261.7 

22000 

432.4 

63.4 

2361.2 

1540.2 

225.6 

8403.1 

22200 

419.2 

48.8 

2216.7 

1583.0 

183.7 

8368.8 

22400 

402.3 

37.2 

2076.4 

1600.5 

148.1 

8261.9 

22600 

378.6 

28.4 

1918.5 

1583.6 

118.4 

8020.1 

22800 

345.1 

21.1 

1723.2 

1499.8 

92.2 

7487.0 
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Illuminant  A 

Illuminant  C 

Eax 

EaV 

Eaz 

Ecx 

Ecy 

Eci 

23000  I  302.0 

15.5  ! 

\ 

1491.2  1 

1352.6  . 

69.2  , 

6677.5 

23200  253.8 

11.1 

1241.6  ! 

1165.1 

49.8  : 

5699.9 

23400  ;  200.6 

7.5 

973.7 

936.9 

34.4  1 

4550.6 

23600  i  142.5 

4.8  ' 

687.5 

669.9 

22.3  ; 

3237.9 

23800  i  94.2 

2.9 

452.0 

445.2 

13.4  ! 

2139.8 

24000  1  62.0 

1.7  1 

295.9 

294.2 

7.9  i 

1406.9 

24200  39.4 

1.1 

188.1 

185.4 

4.9 

885.3 

24400  23.2 

0.7 

111.1 

107.4 

2.9  : 

512.4 

24600  13.6 

0.4 

65.2  1 

62.0 

1.6  1 

294.5 

24800  9.1 

0.4  ' 

44.2  1 

41.3 

1.0  1 

1 

196.4 

25000  6.1  1 

0.3  ! 

29.8  j 

27.2  j 

0.6  1 

129.3 

25200  3.6 

0.2 

17.6 

15.5 

0.3 

73.3 

25400  2.4  I 

0.1  1 

;  11.6 

9.9 

0.1 

46.7 

25600  1.7  i 

7.8 

6.5 

30.6 

25800  1.0  ; 

4.5 

3.5 

16.3 

26000  0.6 

3.0 

2.1 

• 

10.2 

26200  0.4 

i  2.2 

1.5 

7.5 

26400  0.2 

1.4 

0.9 

4.3 

26600  0.1 

1  0.7 

0.4 

2.1 

26800 

0.3 

0.2 

0.9 

27000 

1 

1 

0.1 

'  109828.5  : 

i 

99999.6 

35547.4 

1  98041.3 

100000.0 

118100.4 

for  the  light  of  wavenumber  v  be  f?v,  then  the  tristimulus  values  X,  Y,  Z  of  the 
reflected  light  are  given  by 

A"  000 

li.OM 

Y=”z  R^Efy^vim,  000  ,  (11) 

13. OW 

z=^”zR>Eyi.jvim,ooo 

13.000 

The  problem  under  consideration  is  to  determine  the  interference  color  of  a 
stained  lens  surface  coated  with  anti -reflection  film.  The  program  for  computation 
was  planned  to  calculate  the  tristimulus  values  of  reflected  light  for  various  optical 
thicknesses  of  the  "  stain  ”,  nidi,  the  optical  thickness  of  the  anti-reflection  coating, 
fitdt  being  given. 

To  do  this,  at  first,  sin32  and  cos32  for  the  uppermost  layar  were  calculated 
St  every  calculation  point  and  their  values  were  memorized  by  the  machine  and 
used  in  common  in  the  calculation  for  every  value  of  nidi. 

Following  the  method  outlined  above,  the  calculation  was  performed  using 
FACOM-128  automatic  relay  computer.  FACOM-128  is  a  digital  computer  of 
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floating  system,  built  for  universal  purpose.  The  number  of  effective  flgures  is. 
eight.  The  time  needed  for  the  calculation  in  converting  the  table  of  ExXk  into- 
the  table  of  fvXv  in  equal  intervals  of  v,  was  14  min,  amounting  to  42  min  in  total 
for  .£‘viv,  fvVv,  fv^v.  To  prepare  from  them  an  input  tape  necessary  for  the  next 
step,  it  took  4  min.  Colorimetric  calculation  for  a  given  «idi  took  90  min  for  7 
values  of  n^di. 

The  results  of  computation  will  be  published  in  a  separate  paper”. 

4.  Conclusion 

As  it  is  necessary  for  colorimetric  calculation  to  take  many  calculation  points,, 
it  is  too  laborious  to  do  it  by  hands,  whereas  with  an  automatic  computer,  repeated 
calculations  are  very  easily  performed. 

It  is  not  suited  to  a  mechanical  computer  of  slow  speed  to  calculate  separately 
the  values  of  trigonometrical  functions  required  in  the  determination  of  interference 
color.  Hence  the  program  was  designed  so  as  to  calculate  values  of  the  functions 
by  means  of  the  sum-rule  of  trigonometrical  functions. 

It  was  found  that  the  numerical  integration  using  equal  intervals  in  wave- 
number  is  preferable  to  the  integration  using  equal  intervals  in  wavelength.  It 
also  gives  a  possibility  of  saving  the  labor  in  referring  to  the  trigonometrical  tables 
because  numerical  values  in  function  tables  are  generally  given  at  equal  intervals 
of  variables. 

If  it  is  possible  to  make  use  of  an  electronic  computer  of  large  capacity  for 
memory,  the  device  described  in  this  paper  will  become  unnecessary  and  the 
calculation  inclusive  of  the  change  of  refractive  index  with  the  wavelength  will 
be  performed  easily. 

The  author  wishes  to  express  his  thanks  to  members  of  the  Calculation  Center 
of  the  Fuji  Automatic  Computer,^  who  helped  him  in  performing  calculations 
involved  in  this  work. 


7)  K.  Miyake:  Science  of  Light  6  (1957)  85. 
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Abstract 

When  detection  is  made  by  using  reflection  of  light,  detectable  minimum  of 
optical  thickness  of  stains  formed  on  lens  surface  is  determined  to  be  13-20  m/i. 

As  the  angle  of  incidence  of  the  light  increases,  the  accuracy  falls. 

The  interference  color  of  stained  lens  surfaces  coated  with  anti-reflection 
film  is  determined  using  an  automatic  relay  computer. 

The  results  are  plotted  on  the  CIE  chromaticity  diagram  and  the  color 
difference  is  calculated. 

The  higher  the  refractive  index  of  substrate  glass,  the  smaller  the  detectable 
minimum  of  stain  thickness  when  coated  with  anti-reflection  film. 

Amber  coating  has  an  effect  similar  to  using  a  glass  of  low  refractive  index 
and  may  have  a  possibility  to  prevent  the  appearance  of  irregularity  in  the 
interference  color  of  anti-reflection  coatings. 

1.  Introduction 

In  the  course  of  lens  making,  an  anomalous  layer  called  “  stain  ”  appears  at 
times  on  lens  surfaces. 

To  prevent  the  reflection  of  light  at  lens  surfaces,  nearly  all  lenses  of  high 
grade  are  now  coated  with  anti-reflection  Aims.  If  the  reflection-reducing  coating  is 
deposited  by  evaporation  on  the  lens  surface  on  which  a  stain  has  been  formed, 
the  stain  is  easily  detected  as  irregularity  in  the  color  of  the  coating.  , 

In  optical  workshops,  they  examine  lens  surfaces  generally  with  reflected  light 
and  reject  the  lens  on  the  surface  of  which  any  stain  is  detected.  But  slight 
stains  which  are  unnoticed  in  early  stages,  become  clearly  visible  after  the  coating 
wcis  made  which  is  very  annoying  to  lens  makers. 

By  analysing  the  above  phenomenon,  tolerable  stain  thickness  for  avoiding 
the  otherwise  cumbersome  procedure  was  determined. 

2.  The  Method  of  Detection  Using  Reflected  Light 

Ordinary  method  of  detecting  stains,  adopted  nowadays  in  optical  workshops, 
is  to  let  a  light  from  a  wide  source  such  as  the  blue  sky  or  a  frosted  shade  of 
incandescent  electric  lamp  reflect  at  lens  surfaces  and  to  examine  the  surface.  In 
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this  method,  detection  is  by  the  difference  in  reflectivity  between  the  parts  stained 
and  unstained.  * 

Although  the  accuracy  of  detection  obtained  by  this  method  has  already  been 
calculated  by  Doi^‘,  a  similar  calculation  is  repeated  to  make  a  comparison  with  the 
results  obtained  in  the  following  sections  and  to  enable  a  more  detailed  discussion. 
Assume  that  the  stain  has  a  structure  shown  in  Fig.  1. 


Fig.  1.  Structure  of  the  stain 

The  case  of  normal  incidence  of  light  is  considered  at  first.  Let  ru  r,  be  the 
amplitude  reflection  coefficients  at  the  interfaces  stain— glass  and  air— stain,  n,,  n  the 
refractive  indices  of  glass  and  stain,  /?#,  R  the  amplitude  reflection  coefficients  of 
glass  surface  and  stained  surface.  Then, 


(1) 

r,= 

n  +  ng 

*  1  +  n 

i?_  n  + 

(2) 

•  l  +  r,r, 

where  J={ArclX)nd,  and  d  and  X  being  the  thickness  of  stain  and  the  wavelength  of 
light  considered.  Taking  the  contrast  of  lights  reflecting  at  the  glass  surface  and 
at  the  stained  surface  as 


V= 


Ro  R  * 

/?.!*+  /?!*  ’ 


after  some  calculations,  it  is  given  ,a8 


(3) 


_ 4rirt(l-ri*)(l-rt*)sinM/2  ... 

2(ri  +  rjRl  +  rir,)*  -4rir,  {(ri  +  r,)*  +  (1+ sin* J/2  '  ^ 

Expressing  the  least  contrast  perceptible  to  the  naked  eye  as  V«<„,  the  minimum 

detectable  phase  difference  angle  J  is  given  by 

2  2r»r,  (l-ri*)(l-r,*)+ {(r,  +  rO*  +  (l  +  rir*n  l/.u 

Assuming  ti,  r,,  V*<«  to  be  small  compared  with  unity  and  taking  ri=r,,  we  obtain 
as  an  approximate  expression. 


1)  Y.  Doi:  Tech.  Rep.  Japan  Camera  Res.  Ass.  No.  1  (Nov.  20,  1956)  p.  3  (in  Japanese). 
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sin*|=2V«„  (6) 

The  value  of  Vm<m  depends  on  the  conditions  of  observation.  The  value  known 
as  the  Fechner  fraction  in  photometry,  was  adopted  here.  When  the  average 
intensity  I  is  sufficiently  high,  the  minimum  perceptible  difference  in  intensity  dl 
is  given  as  J///  =  l/60.  As  dljl  is  equal  to  2V»<.,  taking  V«<,=1/120,  the  refrac¬ 
tive  index  of  stain  «=1. 47  and  the  refractive  index  of  the  glass  w,=1.74,  1.63, 
1.52,  the  detectable  minimum  of  optical  thickness  nd  of  the  stain  was  calculated  and 
shown  in  Table  1.  The  wavelength  of  light  is  taken  as  the  one  having  the 
maximum  visual  sensitivity,  that  is,  ^=555m/i. 

Table  1. 

Hf  1.74  1.63  1.52  According  to  (6) 

nd  13.0  14.1  21.1  11.5  m/i 

From  this  result,  it  is  concluded  that  the  detectable  minimum  value  nd  of  stain 
on  glass  used  for  optical  purpose  is  about  from  13  to  20  m^  and  the  higher  the 
refractive  index  of  the  glass,  the  smaller  becomes  this  minimum. 


3.  Effect  of  Incident  Angle 

The  expressions  (1).  and  (2)  are  given  for  the  case  of  normal  incidence.  In 
the  case  of  oblique  incidence,  using  the  Fresnel  coefficients  for  p-  and  s-components 
of  light  as  the  amplitude  reflection  coefficients  in  corresponding  cases,  the  expres¬ 
sions  (2)  and  (4)  remain  unchanged  for  each  of  the  two  components  except  that 
J=(4nlX)ndcos^,  ^  being  the  angle  of  refraction  in  the  stain. 

As  the  result,  the  form  of  the  expression  (5)  remains  the  same  and  using  the 
approximate  expression  (6)  for  ri=r„  the  accuracy  of  detecting^  does  not  change 
in  the  first  approximation.  Therefore  rid  becomes  larger  in  proportion  to  the  factor 
sec^. 

To  verify  this  conclusion,  an  observation  was  made  using  the  apparatus  con¬ 
structed  previously  for  detecting  stains*’.  A  disc  of  SK5  glass  was  mounted  on 
a  supporter  S,  shown  in  Fig.  2.  A  half  of  the  surface  of  the  specimen  was  covered 
with  paraffin  and  the  exposed  half  was  etched  by  dipping  it  in  a  IN  HNO|  aqueous 
solution  at  30°C  for  25  seconds. 

The  glass  surface  was  observed  using  a  light  reflected  at  it.  The  angle  of 
incidence  was  changed  by  the  rotation  of  S  accompanied  with  the  corresponding 


2)  K.  Miyake:  Science  of  Light,  (Tokyo)  5  (1956)  130. 
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Fig.  2 
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change  of  the  direction  of  sight.  The  result  is  shown  in  Fig.  3.  The  contrast 
becomes  lower  with  the  increasing  incident  angle. 

A  conclusion  obtained  from  this  result  is  that  the  smaller  the  angle  of  incidence, 
the  smaller  the  detectable  stain  thickness  provided  that  the  light  used  is  sufficiently 
strong  in  intensity  and  that  no  dbturbing  factor  such  as  reflection  at  back  surface 
inconveniences  the  observation. 

4.  Irregularity  of  Color  of  Anti>reflection  Filins 

The  interference  color  of  thin  films  has  been  investigated  by  Kubota*‘  and 
Monch*’.  Kubota  determined  the  chromaticity  of  color  of  films  deposited  on 
glass  surface  to  various  thickness,  using  the  colorimetric  method  adopted  by  the 
CIE  committee  on  colorimetry. 

If  the  stain  is  so  slight  that,  after  being  coated  with  anti-reflection  film,  there 
is  no  (lerceptible  difference  between  the  interference  color  on  the  stain  and  that 
on  the  rest  of  the  surface,  such  stains  can  be  left  out  of  the  question  in  practice. 
The  maximum  optical  thickness  of  permissible  stain  should  then  be  determined, 
for  which  purpose  the  interference  color  was  examined  for  stains  of  known 
thicknesses. 

The  required  calculation  was  performed  on  an  automatic  relay  computer.  To 
facilitate  the  operation,  the  colorimetric  calculation  was  carried  out  using  the 
wavenumber  as  the  variable,  not  the  wavelength  as  usually  adopted.  Details  of 
the  computation  will  be  published  in  a  separate  paper^’. 

The  structure  of  the  stain  coated  with  anti-reflection  film  was  assumed  to 
be  a  double  layer  shown  in  Fig.  4,  the  refractive  index  «#  of  the  substrate  glass 
being  1.74,  1.63,  and  1.52,  the  refractive  indices  of  the  stain  and  the  reflection- 
reducing  film,  «2  1. 47  and  1. 38  respectively. 


Fig.  4 


3)  H.  KuboU:  J.  Opt.  Soc.  Am.  40  (1950)  146. 

4)  M6nch:  Optik  9  (1952)  97. 

5)  K.  Miyake:  Science  of  Light,  (Tokyo)  6  (1957)  77. 
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As  the  light  source,  the  C-illuminant  of  CIE  was  used  except  in  one  case,  in 
which  the  A-illuminant  of  CIE  was  used. 

Table  2. 

(1)  A-illuminant,  138.9  m/t,  »io=1.74 


«,di  (m/i) 

X 

V 

Y 

0 

0.3934 

0.2757 

0.03059 

1 

0.3832 

0.2723 

0.03028 

2 

0.3732 

0.2693 

0.03008 

4 

0.3538 

0.2643 

0.02999 

9 

0.3120 

0.2588 

0.03156 

16 

O.'ZJM 

0.2652 

0.03806 

25 

0.2597 

0.2865 

0.05343 

36 

C-illuminant,  n2</2= 

0.2694 

138.9m/<,  «o=1.74 

0.3159 

0.08187 

(nift) 

X 

y 

Y 

0 

0.2317 

0.1432 

0.03118 

1 

0.2263 

0.1401 

0.03138 

4 

0.2124 

0.1336 

0.03265 

9 

0.1960 

0.1297 

0.03690 

16 

0.1844 

0.1342 

0.04722 

25 

0.1814 

0.1486 

0.06744 

36 

C-illuminant,  n^2= 

0.1871 

:138.9m;f,  no=l-63 

0.1707 

0.10137 

»\di  (nift) 

X 

y 

Y 

0 

0.2696 

0.2268 

0.06796 

4 

0.2585 

0.2179 

0.06886 

9 

0.2469 

0.2098 

•  0.07144 

16 

0.2348 

0.2036 

0.07772 

25 

0.2258 

0.2026 

0.09003 

36 

C-illuminant, 

0.2224 

:138.9m^,  nfl=1.52 

0.2088 

0.11072 

Hidi  {mft) 

X 

y 

Y 

0 

0.295P 

0.2829 

0.1304 

4 

0.2923 

0.2804 

0.1307 

9 

0.2891 

0.2775 

0.1315 

16 

0.2849 

0.2740 

0.1335 

25 

0.2803 

0.2709 

0.1375 

36 

C-illuminant,  n2</2= 

0.2764 

=  115m;<,  Mq=1.74 

0.2692 

0.1441 

«,di  (m/f) 

X 

y 

Y 

0 

0.4601 

0.3924 

0.07953 

4 

0.4499 

0.3698 

0.06825 

9 

0.4242 

0.3305 

0.05634 

IG 

0.3674 

0.2640 

0.04403 

25 

0.2868 

0.1912 

0.03612 

36 

0.2221 

’  0.1532 

0.03887 
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The  optical  thickness,  nidi  of  the  reflection-reducing  film  was  selected  as 
138. 9mj«  and  115  m/«.  The  former  value  corresponds  to  a  coating  of  sensitive 
purple  color  and  the  latter  to  amber  coating.  The  results  are  tabulated  in  Table  2. 

They  are  plotted  in  Fig.  5,  on  the  CIE  chromaticity  diagram.  In  the  figures, 
white  spots  represent  the  results  of  calculation  made  by  Kubota*’  reduced  to  the 
case  under  consideration  and  black  spots  represent  the  results  of  the  author.  The 
numbers  attached  to  white  spots  mean  the  optical  thickness  of  the  reflection-re¬ 
ducing  film  in  mu  and  those  attached  to  black  spots  are  the  optical  thickness 
Mirf,  of  the  stain. 

The  cases,  in  which  nidi=0,  corrspond  to  the  cases  of  unstained  glass  and 
coincide  with  the  Ctises  investigated  by  Kubota.  Inspecting  the  diagrams,  we 
notice  that  the  coincidence  between  the  results  of  this  work  and  those  of  Kubota 
is  considerably  good. 

Though  Kubota  adopted  an  approximate  expression  by  neglecting  the  effect 
of  multiply  reflected  light,  this  result  shows  that  the  approximation  used  by 
Kubota  is  fair. 

With  the  increase  of  n,di,  black  spots  run  parallel  to  the  locus  of  white  spots. 
This  means  that,  as  the  refrative  index  of  the  stain  does  not  differ  much  from 
that  of  anti-reflection  film,  there  is  a  close  similarity  between  the  stain  and  the 
anti-reflection  film  in  the  effect  of  increasing  the  thickness.  But  if  n^di  becomes 
larger,  the  similarity  is  broken  and  discrepancy  increases  gradually. 

The  above  results  expressed  in  chromaticity  was  converted  to  color  difference, 
taking  account  of  the  difference  in  intensity.  To  calculate  the  color  difference, 
the  method  of  expressing  small  color  differences  studied  by  McAdam”  and  modified 
by  Davidson  and  Friede’’,  and  the  uniform  lightness  and  chromaticity  scale 
(ULCS)  of  Adams*’  were  used. 

The  expression  representing  MacAdam’s  results  is 

jF_rgndJ*+2g^2</jdy+g2tdy*  ,  {dVf  1*'*  ,rj. 

L  (2.5)*  (0.15)*  J  ^  ^ 

<iiu  2gi2,  g22  are  given  graphically  as  functions  of  x  and  y  in  MacAdam’s  paper*’, 

dx  and  dy  are  differences  in  x  and  y  of  the  colors  that  are  compared,  dV  is  a 
difference  in  the  Munsell  V-value. 

According  to  Adams’  method,  the  color  difference  dE  is  given  as 

JE  =  {(0. 23 J  V.)* + [  J(  Vx  -  V,)]* + [0. 4(  F.  -  F,)]*}  *'*  ( 8 ) 

where  F,  is  the  Munsell  F-function,  F,  and  F,  being  the  functions  obtained 

6)  D.  L.  MacAdam:  J.  Opt.  Soc.  Am.  33  (1943)  18. 

7)  H.  R.  Davidson  and  E.  Friede:  J.  Opt.  Soc.  Am.  43  (1953)  581. 

8)  E.  Q.  Adams:  J.  Opt.  Soc.  Am.  (1942)  168. 
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by  replacing  Y  in  V-function  with  Xt  and  Z»,  and  d(V’x— V,)  means  a  difference 
in  (V,— V,)  of  colors  to  be  compared. 

The  color  differences  calculated  by  using  these  expressions  are  shown  in  Fig.  6. 


J£ 


(a)  (b) 

Fig.  6 

The  color  differences  calculated  by  these  two  methods  are  of  similar  tendencies: 
with  a  given  thickness  of  the  stain,  the  more  refractive,  the  larger  the  colof 
difference. 

According  to  a  statistical  investigation  by  Davidson  and  Friede”,  the  color 
difference  perceptible  by  the  naked  eye  with  the  probability  of  50%,  is  given  as 
dE=\  for  MacAdam’s  method  and  JE =0.025  for  Adams’  method. 

Adopting  this  criterion  to  the  case  considered  here  and  doubling  these  limits 
by  considering  the  effect  of  the  light  reflected  at  the  back  surface,  it  seems  likely 
that  the  stains  unnoticed  by  ordinary  method 'of  examination  with  reflected  light 
are  apt  to  become  detectable  by  sensitive  color  of  reflection-reducing  films. 

But  in  optical  workshops,  this  has  been  experienced  only  in  highly  refractive 
glass  such  as  SF3  and  not  in  glass  of  low  refractive  index.  To  explain  this  fact 
satisfactorily,  the  perceptible  limit  of  the  color  difference  must  be  chosen  approx- 
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imately  as  i£=20for  MacAdam’s  method  and  J£=0. 1  for  Adams’  method.  But 
it  may  be  too  arbitrary  to  do  so.  The  values  of  the  perceptible  limit  cited  above 
were  determined  under  a  special  condition  and,  strictly  speaking,  they  are  not 
applicable  to  all  cases.  So  we  must  give  up  the  attempt  to  determine  quantita¬ 
tively  the  detectable  minimum  of  the  thickness  of  stain  coated  with  anti-reflection 
film  using  a  threshold  value  borrowed  from  the  results  of  colorimetric  studies.  The 
crosses  marked  on  every  curve  in  Fig.  6  show  the  detectable  minimum  value  n^dx 
of  stain  obtained  by  using  reflected  light. 

With  the  increase  of  refractive  index  of  glass,  n^di  becomes  smaller,  but  the 
difference  in  interference  colors  becomes  rapidly  noticiable  making  the  detection 
by  anti-reflection  coatings  more  sensitive  than  by  the  ordinary  method. 

This  explains  the  reason  why  faint  stains  may  become  plainly  detectable  on 
highly  refractive  glasses  when  coated. 

We  can  also  conclude  that  the  use  of  the  thickness  of  coating  that  gives 
amber  color  by  interference  has  the  same  effect  as  the  use  of  glass  of  low  refrac¬ 
tive  index. 


5.  Conclusion  and  Acknowledgement 

Detectable  minimum  of  optical  thickness  of  stains  formed  on  lens  surface 
determined  by  using  reflection  of  light  becomes  smaller  as  the  substrate  glass 
becomes  more  refractive.  The  detectable  limit  by  the  adopted  method  is  13-20  mfi. 

By  determining  the  interference  color  of  stained  lens  surfaces  coated  with 
reflection-reducing  film,  it  was  confirmed  that  the  higher  the  refractive  index  of 
the  glass,  the  more  rapidly  grows  the  irregularity  of  the  interference  color. 

The  amber  coating  has  the  same  effect  as  using  the  glass  of  low  refractive 
index  and  may  be  effective  in  preventing  the  appearance  of  irregularity  in  the 
color  of  coatings. 
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Abstract 

Some  problems  on  the  estimation  of  the  absolute  intensity  of  OH  emissions 
in  night  air  glow  are  discussed,  and  a  method  which  seems  best  to  express 
the  absolute  intensity  is  proposed.  In  addition,  a  brightness  standard  con¬ 
structed  for  the  calibration  of  photometer  is  described. 

1.  Introduction 

In  the  previous  paper  under  the  same  title,  Ogawa,  Nakamura  and  Hashizume*' 
reported  about  an  equipment  for  photoelectric  measurement  of  OH  emissions  in 
night  air  glow  assigned  as  a  part  of  the  night  sky  research  in  the  1957-8  Inter¬ 
national  Geophysical  Year'programs  in  Japan.  That  report  gives  only  the  descrip¬ 
tion  of  the  equipment,  and  no  mention  was  made  as  to  the  methods  of  calculation 
of  the  absolute  intensity.  It  would  be  most  practical  if  all  report  on  the  intensity 
measurement  of  the  air  glow  are  made  in  absolute  units.  The  estimation  of  the 
the  absolute  intensity  of  OH  emissions  involves,  however,  more  cumbersome 
problems  than  that  of  atomic  lines  in  the  air  glow. 

In  this  paper  a  method  of  estimation  of  the  absolute  intensity  of  OH  bands 
will  be  discussed,  and  a  brightness  standard  which  is  in  use  in  the  observatories 
at  Shodoshima  (Institute  for  Optical  Research)  and  Maruyama  (Tokyo  Astronomical 
Observatory)  will  be  described, 

2.  Characteristic  Problems  in  the  Estimation  of  Absolute 
Intensity  of  OH  Emissions 

Some  cumbersome  problems  are  inevitable  in  the  estimation  of  the  absolute 
intensity  of  OH  emissions  because  of  the  following  reasons: 

1°  The  spectrum  of  OH  emission  in  air  glow  is  not  of  single  lines,  but  has 

1)  M.  Ogawa,  M.  Nakamura  and  A.  Hashizume:  Science  of  Light  S  (1956)  139. 
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a  band  structure  composed  of  a  great  number  of  lines  extending  over  a  wide 
wavelength  range  from  the  visible  to  the  infrared. 

2°  The  main  part  of  OH  emissions  falls  in  the  near  infrared  region,  where 
photoelectric  sensitivity  of  photo  tubes  is  so  low  that  individual  lines  in  the  bands 
can  not  be  measured  by  ordinary  photoelectric  devices.  Merely  the  total  intensity 
of  a  band  or  several  bands  is  measurable. 

3'  The  relative  intensities  of  the  bands  have  not  yet  been  established.  The¬ 
oretical  values  calculated  under  some  assumptions*  *’  have  been  reported.  These 
values  are  most  reliable  at  present*',  but  some  experimental  results  are  con¬ 
tradictory®’  . 

4°  Studies  of  extinction,  scattering  and  absorption  in  air  are  not  yet  sufficient 
in  the  near  infrared  region  to  give  accurate  corrections  to  the  observed  intensity 
of  the  night  sky  emissions. 

5*  The  two  color  photometer*’  which  is  proved  to  be  effective  for  01  >15577 
line  in  eliminating  the  extra-terrestrial  components  (mainly  light  from  stars  and 
the  zodiacal  light)  is  not  available  for  OH  emissions. 

The  intensity  of  band  emissions  would  be  completely  defined  if  the  intensities 
of  all  individual  lines  are  given.  If,  however,  the  the  relative  intensity  of  individual 
lines  is  known  in  a  band  or  in  a  band  system,  an  expression  of  the  total  intensity 
of  a  band  or  of  a  band  system  would  be  sufficient  to  define  the  intensity. 

'  The  relative  intensity  of  lines  in  a  single  band  corresponding  to  a  vibrational 
transition  is  possibly  defined  by  means  of  the  rotational  temperature  and  the 
relative  intensity  of  the  multiplet  components.  The  rotational  temperature  of  OH 
bands  in  night  air  glow  was  measured  by  several  authors'".  Meinel  *  gave  265' K 
as  the  rotational  temperature  and  0. 317  as  the  intensity  ratio  of  *ni/i  and  *IIi/» 
from  his  observations  at  a  moderate  latitude. 

The  deflections  of  photometer  are  due  not  only  to  the  air  glow  emissions  but 
also  to  the  light  from  stars  ai^d  to  the  zodiacal  light.  These  extra-terrestrial 


2)  H.  S.  Heaps  and  G.  Herzberg:  Zeits.  f.  Phys.  133  (1952)  48. 

3)  F.  E.  Roach:  The  International  Geophysical  Year  Instruction  Manual  IV,  Pergamon 
Press,  London  (1956),  p.  76. 

4)  A.  V.  Jones:  Nature  172  (1953)  496. 

H.  P.  Gush  and  A.  V.  Jones:  J.  Atm.  Terr.  Phys.  7  (1955)  285. 

5)  N.  I.  Eedorova:  Astronomchski,  Zhurnal,  34  (1957)  1247. 

6)  F.  E.  Roach  and  H.  B.  Pettit:  J.  Geophys.  Res.  56  (1951)  325. 

M.  Huruhata,  T.  Nakamura  and  H.  Tanabe,  Tokyo  Astronomical  Bulletin,  Sec.  Ser. 
9A  (1957)  1005. 

7)  A.  B.  Meinel:  Astrophys.  J.  112  (1950)  120. 

J.  Cabannes,  J.  Dufay:  Compt.  Rend.  232  (1951)  426. 

C.  Oliver:  Phys.  Rev.  90  (1953)  1118. 
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components  must  be  substracted  from  the  original  values  of  deflections  of  the 
photometer.  If  an  assumption  that  the  spectra  of  the  star  light  and  of  the  zodiacal 
light  are  both  Go  type  in  average  is  allowed,  one  can  estimate  the  intensity  of 
the  extra-terrestrial  components  at  a  wavelength  in  the  infrared  region  from  those 
at  5300  A,  at  which  the  intensities  are  easily  measured  because  of  the  absence  of 
contribution  of  the  air  glow.  The  variation  of  the  extra-terrestrial  components 
with  the  zenith  distance  is  known  to  be  small  at  5300  A  which  would  also  be  the 
case  with  the  infrared  region.  The  effect  of  absorption  by  water  vapor  in  the 
extra-terrestrial  components  calculated  from  the  amount  of  precipitable  water  and 
the  absorption  coefficients  of  water  vapor  is  only  6%  in  summer  and  3%  in  winter 
for  our  photometer.  As  the  actual  contribution  of  the  extra-terrestrial  components 
in  photometer  deflection  in  the  near  infrared  region  is  only  about  20%,  somewhat 
rough  estimation  of  this  background  would  be  tolerable. 

From  above  reasons,  it  would  be  possible  to  give  the  proper  deduction  of  the 
extra-terrestrial  components  from  the  deflections  of  the  photometer. 

The  effects  of  extinction,  scattering  and  absorption  by  the  earth’s  atmosphere 
for  the  air  glow  components  can  be  eliminated  by  means  of  ”  the  mean  value 
method  ”  in  which  the  multiplication  factors  for  observed  values  differing  in  zenith 
distance  are  found  from  their  mean  values.  This  is  based  on  the  reasoning  that 
the  mean  of  air  glow  intensities  observed  at  one  place  is  the  same  at  any  other 
place,  provided  that  ob^rvations  are  made  for  a  long  period. 


3.  Calculation  of  the  Abaolnte  Intensity  of  OH  Emissions 
The  deflection  of  the  photometer  will  be  given  by  the  following  formula. 

D = A  = A  J^G(>l)5(>l)d/l , 

where  D:  deflection  of  the  photometer, 

A:  sensitivity  of  electronic  accessory  of  the  photometer, 

S(il):  spectral  sensitivity  of  the  photo  tube  used, 

Tfi):  overall  spectral  transmission  of  the  optical  elements, 

G(>i)=SW-T(A), 

B{i):  brightness  of  the  night  sky, 

X:  wavelength. 

B(i)  can  be  divided  into  two  components,  namely  the  air  glow  emissions  and 
the  extra-terrestrial  lights. 

B{i)=B^M.O)+B,.t.O) 


As  the  air  glow  emissions  consist  of  a  number  of  lines. 
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D=A  + i4  t.  , 

=  Dm.K.  +  D,.t.  , 

where  Bm.t.W  means  the  integrated  brightness  over  a  spectral  line  or  a  narrow 
spectral  band.  Both  SW  and  T{X)  and  therefore  G{i)  can  be  determined  from 
laboratory  measurements.  .Be.t.(>l)  can  be  calculated  from  photometric  data  at  5300  A 
and  the  spectral  distribution  of  the  extra-terrestrial  coponents.  A  can  be  determind 
by  using  a  brightness  standard  of  known  spectral  distribution  and  absolute  intensity. 
If  the  standard  emits  only  a  monochromatic  or  a  nearly  monochromatic  radiation 
of  wavelenght  the  deflection  of  the  photometer  for  this  standard  is 

D^— AG{X—Xi^I}^  ^  ' 

or  A— DJG{X=Xt)Bt  ^ 

where  B^  is  the  known  brightness  of  the  standard.  The  extra-terrestrial  part  in 
the  photometer  deflections  is,  therefore, 

As  the  relative  intensities  of  spectral  lines  or  bands  in  air  glow  would  not  vary 
with  time  or  place,  we  can  put 

Bu.tXX)=E-FW, 

where  E:  a  constant, 

F{X):  relative  intensity  of  the  spectral  components  (lines  or  bands)  in  the 
air  glow. 

Then, 


or 


At  present,  relative  intensities  F(i)  of  OH  bands  are  not  known  for  certain. 
It  seems  best  to  express  the  absolute  intensity  of  OH  emission  in  the  air  glow 
provisionally  by  (D— De.t.)Bo/^o  and  G(/l),  since  Bm.uXX)  values  can  be  easily  cal¬ 
culated  from  the  above  equation  after  the  exact  F{i)  values  are  found. 


4.  Brightness  Standard 


A  brightness  standard  was  constructed  to  determine  the  absolute  sensitivity 
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of  the  photometer.  This  standard,  which  is 
shown  in  Fig.  1,‘can  be  easily  mounted  on  the 
air  glow  photometer  by  means  of  a  shoe.  The 
light  source  of  the  standard  is  a  midget 
tungsten  lamp  which  is  lighted  at  a  current 
of  110  ma  dc,  while  its  normal  rate  is  300  ma. 

The  lamp  illuminates  a  disk  of  opal  glass 
which  is  placed  40  cm  apart  from  the  lamp. 

An  interference  filter  of  maximum  transmission 
at  8900  A  was  inserted  between  the  lamp  and ' 
the  glass  disk.  The  half  width  of  the  spectral 
transmission  of  the  filter  is  300  A.  The  opal 
glass  disk  is  12  cm  in  diameter  and  is  large 
enough  to  cover  the  entire  field  of  view  of  the 
photometer. 

The  brightness  of  the  opal  glass  was 
calibrated  by  means  of  a  thermocouple  and  a 
Mazuda  radiation  standard  lamp.  The  bright¬ 
ness  of  the  standard  by  110  ma  current  is 
2. 61 X 10"*  microwatts/cm*  steradian,  which  corresponds  to  1. 17  x  10**  quanta/sec 
cm*  steradian  for  photons  of  8900  A,  or  to  1. 47  x  10®  rayleighs  for  the  air  glow  * 
intensity. 

5.  Absolute  Intensity  of  OH  Emissions  in  the  Air  Glow 

The  intensity  of  OH  emissions  in  the  air  glow  greatly  varies  with  time  and 
direction  viewed.  The  deflection  of  the  photometer  is  from  about  one  tenth  to  a 
quarter  of  that  for  the  standard,  when  it  is  directed  to  zenith  in  ordinary  clear 
nights.  The  extra-terrestrial  components  are  estimated  to  be  about  one  thirtieth 
of  the  brightness  standard,  so  the  deflection  due  to  the  air  glow  alone  is  from 
one  fifteenth  to  one  fifth  of  that  of  the  standard. 

Then 

D  1 

15~  5  ’ 

So =1.47x10®  rayleighs  at  8900  A, 

DBJDt = (1. 0~3. 0)  X 10®  rayleighs. 

If  the  relative  intensities  of  OH  bands  given  by  Roach  in  the  IGY  instruction 
manual*’  are  used. 
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G(^=8900A)=5.0, 

F(7-3  band) =47, 

F(>l)GW)=1998. 

The  absolute  intensity  of  7-3  band  of  OH  emissions  of  the  air  glow  is,  therefore, 
(1. 2~3. 6)  X 10*  rayleighs, 

which  agrees  with  Roach’s  value  of  1940  rayleighs*'. 

6.  Summary 

Some  problems  on  the  estimation  of  the  absolute  intensity  of  OH  emissions 
in  night  air  glow  were  discussed,  and  a  method  of  expressing  the  absolute  intensity 
was  proposed.  In  addition,  a  brightness  standard  constructed  for  the  calibration 
of  photometer  was  described. 
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Abstract 

The  influence  of  addition  of  flux  on  the  luminescent  properties  and  ultra¬ 
violet  absorption  of  Cu-activated  CaS  phosphors  is  investigated.  As  flux, 
sulfates,  fluoride  and  chlorides  of  Li,  Na  and  K  are  used.  According  to  the 
results  of  experiments,  Cu-activated  CaS  phosphors  show,  in  general,  two 
emission  bands  at  about  430  m/i  and  495  m/i  and  the  sulfate  and  fluoride  fluxes 
intensify  the  latter  band,  whereas  chloride  fluxes  the  former.  In  addition,  it 
may  be  inferred  that  the  latter  band  originates  from  CaS  host  crystal  itself 
and  the  former  from  Cu-activator  located  in  it. 

1.  Introduction 

The  effect  of  the  flux  such  as  NaCl  and  KCl  in  Cu-activated  ZnS  phosphors 
has  been,  to  some  extent,  clarified  by  several  investigators'*  According  to  the 
recent  investigations  the  pair  of  Cu-ion  and  Cl-ion,  which  are  located  substitu- 
tionally  in  CdS  host  crystals  plays  the  essential  role  also  for  the  photoconductivity 
of  CdS  crystals'*’.  On  the  other  hand,  for  CaS  phosphors  sulfates  and  fluorides  of 
alkali-metals  such  as  Na2S04,  LiF,  NaF  etc.  have  been  used  customarily  as  flux*'. 
However,  effect  of  these  fluxes  on  the  Cu-activated  CaS  phosphors  has  not  been 
completely  clarified.  This  paper  deals  with  the  influence  of  flux  added  to  Cu- 
activated  CaS  phosphors  on  their  visible  luminescence  and  ultraviolet  absorption. 

2.  Structure  and  Composition  of  host  crystal.  ' 

The  phosphors  submitted  to  the  experiments  were  prepared  by  Leriard’s 
method,  i.e.  they  are  synthesized  by  heating  the  mixture  of  pure  sulfur  and  CaO 
powder  at  950“C  for  20  min.  in  an  electric  furnace  after  the  addition  of  adequate 
quantities  of  flux  and  activator.  As  the  activator  CUSO4  and  as  the  flux  sulfate, 

*  The  results  presented  in  this  jsaper  were  briefly  reported  at  the  Autumn  Meeting,  held 
.at  Okayama  in  October  1957,  of  the  Society  of  Applied  Physics  (Japan). 

1)  F.  A.  Kroger  and  J.  E.  Hellingman:  J.  Electrochem.  Soc.  93  (1948)  156,  95  (1949)  68. 

2)  R.  H.  Bube:  J.  Chem.  Phys.  19  (1951)  985. 

3)  R.  H.  Bube:  J.  Chem.  Phys.  20  (1952)  708. 

4)  F.  A.  Krdger  et  al:  Physica  15  (1949)  990. 

5)  R.  H.  Bube  and  S.  M.  Thomsen:  J.  Chem.  Phys.  23  (1954)  15,  18. 

6)  F.>  Schmidt:  Ann.  d.  Phys.  19  (1932)  211. 
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chloride  or  fluoride  of  Li,  Na  or  K  were  used.  Proportions  of  activator  and  flux 
added  are  shown  in  weight  percent. 

Sulfur  was  purified  by  thrice-repeated  sublimation  and  purified  CaO  powder 
was  prepared  from  CaCOj’"*’.  In  heating  the  mixture,  double  crucible  was  used 
to  avoid  the  oxidation. 

According  to  the  current  view,  in  Lenard’s  method,  host  crystals  are  composed 
during  the  heating  by  following  reaction*’. 

4S-l-4CaO - ►XaS+CaSO.*  (2.1) 

In  this  method  fairly  large  quantity  of  CaS04  is  inevitably  contained  in  the  host 
crystals  in  addition  to  CaS.  X-ray  inspection  by  powder  method  showed  that  the  host 
crystals  are  the  mixture  of  regular  CaS  and  CaS04  crystals,  regardless  of  the  sort 
and  quantity  of  added  flux.  Further,  any  trace  of  CaO  crystal  was  not  recognized 
with  the  sample  prepared  under  the  above  mentioned  conditions.  In  order  to 
determine  the  ratio  of  the  amount  of  CaSOt  to  the  total  quantity  of  the  host 
crystals,  the  amount  of  sulfate  radicals  in  the  prepared  phosphors  was  measured 
by  chemical  analysis.  5g  of  the  phosphor  is  suspended  in  20  cc  of  hot  acetic  acid 
(cone,  acetic  acid  1 :  water  4).  In  this  procedure  CaS  crystals  in  the  mixed  crystal 
are  decomposed  by  forming  H,S  gas,  whereas  most  part  of  CaS04  crystals  remain 
suspended  in  liquid.  This  turbid  liquid  is  introduced  into  500  cc  of  boiling  water. 
After  stirring  for  a  long  time,  the  CaS04  crystals  may  be  dissolved  gradually  and 
clear  liquid  is  obtained. 

The  amount  of  sulfate  radicals  in  the  clear  liquid  may  be  determined  by  mak¬ 
ing  the  sulfate  radicals  precipitate  as  BaS04*®’. 

The  results  obtained  are  shown  in  Table  I.  In  Table  I  the  quantity  of  sulfate 
radicals,  which  originate  from  activator  and  flux  (when  sulfate  of  alkali-metal  is 
used  as  flux)  are  subtracted  from  the  results  of  analysis.  Referring  to  the  same 
table,  we  can  point  out  that  the  proportion  of  sulphate  radical  in  the  host  crystals 
is  in  good  agreement  with  the  value  26.5%  expected  from  the  reaction  (2.1),  re¬ 
gardless  of  the  sort  and  quantity  of  the  flux  used. 

From  the  results  stated  above,  we  can  infer  that  the  optical  properties  of  pure 
CaS  phosphor  may  be  known,  if  from  the  optical  properties  of  CaS  phosphor  pre- 

7)  P.  Lenard  and  V.  Klatt:  Wied.  Ann.  38  (1899)  90,  Ann.  Phys.  15  (1904)  225,  425,  633. 

8)  P.  Lenard.  F.  Schmidt  and  R.  Tomaschek:  Phosphoreszem  und  Fluoreszenz,  Hd.  d. 

Exp.  Phys.  Vol  23  Part  1  (Akad.  Verlag,  Leipzig  1928)  p.  327. 

9)  P.  Lenard:  Reference  8)  p.  325. 

In  practice,  about  twice  of  the  sulphur  amount  required  to  perform  the  reaction  (2.1) 
are  mixed  with  CaO  powder.  The  excess  of  sulphur  is  sublimated  during  heating. 

10)  W.  W.  Scott:  Standard  Methods  of  Chemical  Analysis  (D.  Van  Nostand  Co.  N.  Y. 

1939)  Vol.  1,  p.  127. 
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Table  I.  Amount  of  S04^~-radical  in  host  crystal. 


Flux 

Amount  of  flux 
(%) 

Amount  of 
S04*--radical  (%) 

Na*S04 

6.00 

27.2 

Na2S04 

0.6 

27.6 

Na2S04 

0.06 

26.8  j 

K2SO4 

0.5 

26.4  j 

KF 

2.0 

27.5 

LiF 

3.2 

25.2 

NaCl 

2.0 

25.4 

pared  by  Lenard’s  methcxi  those  of  CaS04  crystal  prepared  under  the  same  con¬ 
ditions  are  substracted,  so  far  as  the  optical  properties  of  CaS  and  CaS04  crystals 
in  the  host  crystals  are  independent  of  one  another.  For  this  purpose,  CaS04 


powder  crystals  containing  the 
same  quantities  of  activator 
and  flux  as  the  individually 
prepared  CaS  phosphor  were 
prepared. 

3.  Emission  bands. 

Pure  host  crystal  consisting 
of  CaS  and  CaS04  is  nearly 


non-luminescent  under  2537  A 
excitation.  The  phosphor  con¬ 
taining  Na2S04  as  flux  and  no 
copper  shows  an  emission  band 
at  about  495  m^i  which  in¬ 
creases  in  radiance  with  increas¬ 
ing  proportion  of  added  flux  as 
shown  in  Fig.  1  a.  The  phos¬ 
phor  containing  copper  and 
no  flux  shown  two  weak  emis¬ 
sion  bands  at  about  430  m/i 
and  495  m/i.  When  both  flux 
and  copper  are  added  to  the 
host  crystal,  the  two  emission 
bands  are  also  observed.  In 
general,  the  band  at  430  m/i 
increases  in  radiance  more 


400  500  600 


Wave-length  (m/<) 

Fig.  1  a  and  b.  Spectral-distribution  curves  of  fluo¬ 
rescence  of  CaS  phosphors,  excited  by  2537  A  line. 

(1)  Na2S04  0.064% 

(2)  Na8S04  2.000% 

(3)  Cu  6xlO-*%-l-NaCl  0.06% 

(4)  Cu  6x10  *% -HNaCl  0.6% 

(5)  Cu  6xlO-*%-l-NaCl  3.2% 

(6)  Cu  6xl0~»%-^Na2SO4  0.064% 

(7)  Cu  6  X 10-*%  -I-  Na2S04  0 . 64% 

(8)  Cu  6x10  »%-»-Na2S04  3.2% 
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rapidly  than  the  band  at  495  with  increasing  proportion  of  flux,  when  the  flux 
added  is  a  chloride  (cf.  Fig.  la).  On  the  contrary,  the  latter  band  increases  in 
radiance  more  rapidly  than  the  former  with  increasing  proportion  of  flux,  when 
a  sulfate  or  fluoride  is  added  as  flux  (cf.  Fig.  lb).  The  ratio  of  peak  intensity 
/430  of  the  band  at  430  m/i  to  that  Itts  of  the  band  at  495  m/i  under  the  same  ex¬ 
citation  condition  is  plotted  in  Fig.  2  against  the  quantity  of  added  flux. 

KC1(1.2) 

.NaCI(2) 

LiCI  (^5) 


K,SO.(8) 

LiK( 130) 

KF(16) 

NaK(160) 

NajSO.dOO) 

Li.SO.(120) 


005  0.1  05  1  5 

Proportion  of  flux  (%) 

Fig.  2.  Ratio  of  peak-intensities  of  two  emission  bands  versus  proportion  of  added  flux 
under  2537  A  excitation.  All  the  phosphors  contain  6  X 10“*%  of  Cu. 

The  numerical  values  in  the  brackets  appended  to  the  names  of  fluxes  show 
the  relative  magnitude  of  of  the  samples  containing  5  weight  percents  of  flux. 
Referring  to  these  numerical  values,  we  may  point  out  that  Na*-  and  Li^-ions  are 
more  effective  to  luminescence  than  K^-ion,  so  far  as  the  phosphors  with  the 
fluxes  containing  the  same  anion  dr  anionic  radical  are  concerned.  Among  the 
phosphors  containing  the  same  alkali  cation  as  flux,  those  containing  fluoride  and 
sulfate  are  more  efficiently  luminescent  than  those  containing  chloride. 

All  CaSO*  powder  crystals  containing  copper  <md  one  of  these  fluxes  were 
non-luminescent  or  showed  only  negligible  luminescence.  Therefore,  it  may  be 
inferred  that  the  properties  stated  above  on  the  luminescence  may  be  attributed 
to  CaS  crystals  mixed  in  the  host  crystal  and  CaS04  crystals  coexisting  with  them 
have  no  contribution  to  luminescence. 

4.  Ultraviolet  Reflectance 

Ultraviolet  reflection  spectra  of  some  of  the  samples  prepared  are  shown  in 
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Fig.  3.  Pure  MgO  powder  is  used  as  the  standard  of  reflectance  (reflectance: 

The  reflectance  curves  of  all  the  prepared  CaS04  powders  are  nearly 
flat  over  the  wave-length  range  240~4(X)m/x  (Curve  (1)  and  (2)  in  Fig.  3). 


250  300  3S0  400 

Wav«- length  (m/<) 

Fig.  3.  Spectral-reflectance  curves  of  CaS04  powder  crystals  and 
some  CaS  phosphors,  relative  to  MgO  powder. 

(1)  pure  CaS04 

(2)  CaS04  (Cu  6xl0-»?^  +  Na2SO4  2%) 

(3)  CaS  phosphor  (host-crystal) 

(4)  -  n  (Cu  8x10-*%) 

(5)  //  (Cu  6xlO-»%-t-LiF  0.6%) 

I  (6)  //  (Cu  6xlO-»%-fLiF  3.0%) 

Therefore,  we  may  infer  that  CaS04  crystals  in  the  host  crystal  have  no  con¬ 
cern  with  the  wave-length  characteristics  of  the  reflectance  curve  of  the  host 
crystal.  Referring  to  curve  (3)  in  Fig.  3,  we  find  that  the  phosphors  containing 
no  copper  and  flux  shows  a  weak  absorption  band  at  about  260  m^.  When  copper 
is  added  as  activator,  a  strong  absorption  band  appears  near  the  band  stated  above 
(curve  (4)  in  Fig.  3).  The  phosphors  containing  both  copper  and  flux  show  broad 
and  nearly  uniform  absorption  in  the  wave-length  range  3(X)~4(X)m^,  regardless 
of  the  sort  of  the  flux  added,  in  addition  to  the  strong  band  at  260  m/i.  This  uni¬ 
form  absorption  gradually  increases  with  increasing  proportion  of  flux,  whereas 
the  band  at  260  m/u  remains  constant  (curve  (5)  (6)  in  Fig.  3). 

5.  Decay  of  Phosphorescence 

The  decay  of  phosphorescence  was  measured  on  some  of  the  prepared  phos¬ 
phors  by  using  an  electromagnetic  oscillograph  for  early  stage  of  the  decay  and 

11)  .Jh.  P.  J.  Botden  and  F.  A.  Kroger:  Physica  15  (1949)  749. 

12)  S.  Asano:  Sci.  of  Light  (Tokyo)  5  (1956)  66. 
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a  sensitive  galvanometer  for  later 
stage  of  the  decay.  Fig.  4  illustrates 
the  decay  curves  of  the  prepared 
samples  in  log-log  plotting.  Refer¬ 
ring  to  the  same  figure,  we  find  that 
the  slope  of  the  curve  approaches 
to  a  constant  value  for  large  decay 
time.  Then  the  intensity  of  the 
phosphorescence  B{t)  may  be  written 
in  the  form 

B(/)=5o/(l+a/)-,  (5.1) 

where  /  is  the  decay  time  after  the 
cessation  of  excitation  and  n,  a  and 


Decay  time  (lec) 

Fig.  4.  E>ecay  curves  of  green  emission  band 
after  2537  A  excitation  for  3  minutes. 

(1)  Cu  6xlO-*%+UF  3.096 

(2)  Cu  6xl0-»9^-(-Na2SO4  3.2% 


Bo  are  constants  independent  of  t. 
The  values  of  these  constants  should 
be  proper  to  every  particular  sample, 
so  far  as  the  excitation  conditions 


Decay  tim'-  (tec) 

Fig.  5.  Decay  curuves  of  phosphorescence  in 
early  stage. 


i.e.  kind  and  intensity  of  excitant, 
duration  of  excitation  and  operating 
temperature  are  fixed  constants.  The 
values  of  exponent  n  obtained  from 
the  straight  parts  of  the  curves  in 
Fig.  4  are  also  shown  in  the  same 
figure.  Furthermore,  the  aspects  of 
the  earlier  stage  of  these  decay 
curves  are  shown  in  Fig.  5  by  semi¬ 
log  plotting.  The  curves  show  con¬ 
spicuous  bends  about  0. 01  sec.  after 
the  cessation  of  UV  excitation.  The 
shapes  of  the  rapidly  decaying  i>arts 
before  the  bends  are  indefinite,  since 
quick  following  of  the  vibrator 
within  0.01  sec.  was  prevented  by 
its  proper  frequency  (150  cycles/sec.). 
However,  we  may  determine  fairly 


accurately  the  values  of  a  and  Bo  from  the  more  slowly  decaying  parts  of  the 


decay  curves  in  Fig.  5  and  from  the  values  of  n  already  known.  The  values  of 


a  thus  obtained  are  shown  in  the  same  figure. 
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If  we  extrapolate  the  slowly  decaying  parts  in  Fig.  5  towards  shorter  decay 
times  by  using  (5. 1)  and  the  numerical  values  of  n,  a  and  fio  which  are  determined 
in  the  above  manner,  we  obtain  the  dotted  curves  shown  in  the  same  figure. 
These  curves  show  that  the  phosphorescence  consists  of  a  very  quickly  decaying 
one  and  a  long  persistent  one.  It  is  assured  by  preliminary  spectroscopic  observa¬ 
tion  that  the  latter  part  corresponds  to  the  blue  band  at  495  m/<  zmd  that  the 
phosphorescence  of  the  blue  band  at  430  m^*  cannot  be  observed  visually.  There¬ 
fore,  the  former  quickly  decaying  part  will  be  able  to  be  identified  with  the  phos¬ 
phorescence  of  the  blue  band,  although  the  type  of  the  decay  curve  could  not  be 
•determined  in  the  present  experiments. 

6.  Summary 

From  the  results  obtained  in  the  preceding  sections,  the  following  facts  and 
inferences  may  be  pointed  ont. 

(1)  The  host  crystals  of  CaS  phosphor  prepared  by  Lenard’s  method  are  the 
mixture  of  CaS  and  CaS04  powder  crystals  with  the  stoichiometric  proportion 
determined  from  the  reaction  (2.1). 

(2)  Among  the  main  emission  bands  at  495  m/i  and  430  m/4,  the  former  green 
band  will  be  attributed  to  the  host  crystal  CaS  itself,  whereas  the  latter  blue  band 
to  the  copper  activators  in  CaS  crystal,  (cf.  Section  3) 

Furthermore,  the  uniform  absorption  in  the  range  300~400  mti  is  related  with 
the  green  band  and  absorption  at  about  260  m/i  with  the  blue  band.  (cf.  Section  4) 

(3)  The  sulfate  and  fluoride  fluxes  intensify  mainly  the  green  band,  whereas 
the  chloride  fluxes  the  blue  band.  It  should  be  noted  that  Cl'-ion  is  effective  to 
the  intensification  of  the  emission  band  of  copper  in  CaS  crystal  in  a  similar 
manner  as  in  the  case  of  Cu-activated  ZnS  phosphors^'"^’. 

(4)  In  general,  compounds  of  F  and  Na  are  more  effective  to  intensificatidn 
of  luminescence  than  those  of  K.  It  would  be  supposed  that  K^-ion  has  too  large 
ionic  radius  to  occupy  a  site  effective  to  intensification  in  CaS  host  crystal. 


SCIENCE  OF  LIGHT 


VOLUME  6.  NUMBER  3.  1957 


Solvent  Effect  on  O-H  Stretching  Vibration  of  Methanol.  I. 

Yoshiki  Sato 

Institute  for  Optical  Research,  Tokyo  University  of  Education 
Shinjuku-ku,  Tokyo 

(Received  Nov.  15,  1957) 

Abstract 

In  carbon  tetrachloride  solutions  of  various  oxygenated  and  nitrogenated 
organic  compounds,  detailed  measurements  of  the  first  overtone  band  of 
-  methanol  hydroxyl  stretching  vibration  are  made  with  an  infrared  grating 
spectrometer,  which  resulted  in  obtaining  two  new  absorption  bands  in  longer 
wavelengths  of  monomolecular  OH  band.  Predominant  behaviors  are  found 
between  their  intensity  ratios  and  the  basicity  constants  of  the  proton  accepting 
molecules  and  also  between  their  separations  and  the  latter,  which  are  explained 
on  the  basis  of  a  potential  field  consisting  of  two  minima  derived  in  the  previous 
paper  for  the  motion  of  bonded  hydrogen  atom.  Furthermore  the  feature  of 
the  potential  curve  is  likely  made  clear  by  analyzing  the  above-mentioned 
experimental  results. 

1.  Introduction 

In  the  previous  paper*’,  it  was  shown  that  in  infrared  studies  on  the  associa* 
tion  between  substituted  phenols  and  some  proton  acceptors  through  hydrogen 
bond,  two  new  absorption  bands,  both  likely  to  be  ztscribed  to  the  hydrogen 
bonding,  are  observed  in  1. 5ju  region,  the  appearance  of  which  can  be  explained 
by  assuming  a  double  minimum  potential  for  the  motion  of  ixroton.  Recently, 
measurements  in  1.1  region  for  some  of  the  above-described  hydrogen  bonding 
systems  also  resulted  in  obtaining  two  new  absorption  bands,  both  of  which  were 
attributed  to  the  hydrgen  bonding  by  considering  how  the  potential  curve  for  the 
OH  valence  vibration  are  deformed  in  the  complexes  formed  by  the  solvent  mol¬ 
ecules  with  the  solute  molecule.  Further,  in  a  few  lower  levels,  analogous  be¬ 
haviors  of  their  wave  number  differences  to  the  magnitudes  of  the  doublet  splitt¬ 
ings  of  vibration  of  ammonia  allowed  us  to  conclude  that  the  height  of  the 
potential  hill  for  proton  transfer  is  probably  about  20  Kcal./mole,  for  example,  in 
hydrogen  bonding  between  phenol  and  acetone^’. 

In  order  to  throw  further  light  upon  the  above-mentioned  idea,  the  examin¬ 
ing  of  the  dependences  of  intensities  and  separations  of  such  two  absorption  band& 
on  the  basic  properties  of  various  proton  acceptor  molecules  appears  to  be  of 

1)  Y.  Sato  and  S.  Nagakura:  Science  of  Light  [Tokyo]  4  (1955)  120. 

2)  Y.  Sato:  Science  of  Light  [Tokyo]  6  (1957)  69. 
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great  interest,  from  which  some  informations  as  to  the  feature  of  the  potential 
curve  under  consideration  may  be  obtained.  For  thb  purpose,  studies  on  the 
association  through  the  hydrogen  bond  between  methanol  and  several  proton  ac¬ 
ceptors  were  made  in  detail  by  measuring  the  first  overtone  band  of  methanol 
hydroxyl  stretching  vibration  in  carbon  tetrachloride  solutions.  As  a  result,  close 
connections  between  the  quantities  in  question  were  obtained,  which  led  us  to 
determine  approximately  the  feature  of  the  double  minmum  potential  curve  in  which 
the  proton  of  the  associated  donor  molecule  is  in  motion. 

2.  Experimental  and  Results 

Methanol  was  used  as  proton  donor  because  of  its  weak  acidity  and  the  well 
known  0-H  vibrational  frequencies  in  gaseous  state,  and  various  oxygenated  zmd 
nitrogenated  compounds  as  proton  acceptor.  All  the  substances  used  are  of  special 
grade  available,  some  of  which  were  prepared  by  redistillations  before  being  used. 
A  near  infrared  grating  spectrometer,  previously  described  in  earlier  papers,  was 
maintained  at  approximately  constant  conditions  throughout  this  experimental 
work,  with  which  all  spectral  measurements  were  made  in  carbon  tetrachloride 
solutions  of  the  proton  acceptor  molecules  by  recording  spectral  curves  of  methanol, 
from  each  of  which,  as  the  background,  the  absorption  band  of  methanol  itself 
was  substracted.  All  of  values,  obtained  in  this  way,  were  averaged  over  several 
runs. 


Fig.  1.  Absorption  curves  of  methanal  in  carbon  tetrachloride  solutions  of:  (1) 
1.92  mole/1,  acetophenone;  (2)  2.07  mole/1,  diethyl  ether;  (3)  1.37  mole/1,  a-picoline; 
(4)  0.78  mole/1,  tributylamine.  The  dotted  curve  is  for  0.16  mole/1,  methanol  itself. 


Some  of  the  absorption  curves  obtained  are  shown  in  Fig.  1,  in  which  two 
new  absorption  bands  were  found  in  1. 5  region,  their  shifts  from  the  unassociated 


Fig.  2.  Correlations  of  intensity  ratios  and  separations  of  the  two 
bands  with  the  basicity  constants  of  proton  acceptors. 


OH  band  of  methanol  being  listed  in  Table  1.  As  clearly  seen  in  Fig.  1  and 
Table  1,  the  relative  intensities  and  the  separations  of  the  two  absorption  bands 
vary  characteristically  with  the  changing  of  acceptor  molecules.  Then  plotting 
the  ratios  of  the  maximum  optical  densities  or  extinction  coefficients  of  the  longer 
wavelength  bands  to  those  of  the  shorter  ones  and  their  separations  against 
logarithms  of  the  basicity  constants  of  the  acceptor  molecules.  Fig.  2  is  obtained. 

J£ 
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Table  1.  Relations  of  proton  acceptors  to  wave  numbers,  separations  of  the  two 
absorption  bands,  separations  of  the  unperturbed  levels,  their  shifts  from  the  unas¬ 
sociated  OH  band,  and  coefficients  in  mixtures  of  zero  approximation  eigenfunctions. 


Accepter 

(cm"^) 

(cm"') 

JE 

(cm"') 

J  (cm"')< 

Wu 

(cm"') 

(cm"') 

(cm"') 

Acetophenone 

6453 

7004 

551 

0.24 

338 

218 

227 

567 

Diethyl  ketone 

6443 

6968 

525 

Ethyl  acetate 

6435 

6951 

516 

Dioxane 

6431 

6920 

489 

0.47 

176 

228 

361 

537 

Diethyl  ether 

6421 

6859 

438 

0.72 

71 

216 

449 

520 

Quinoline 

6371 

1  6768 

411 

1.69 

102 

192 

606 

504 

a-Picoline 

6357 

6768 

397 

1.75 

112 

198 

618 

506 

Pyridine 

6343 

6747 

405 

1.84 

120 

193 

1  640 

520 

Tributylamine 

6279 

j  6757 

478 

2.46 

202 

217 

708 

506 

Diethylamine 

6240 

6772 

532 

2.54 

231 

240 

734 
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The  basicity  constants  used  here  are  the  dissociation  constants  of  the  bases  for 
water  solution  divided  by  the  autoprotolysis  constant  of  water,  which  were  taken 
mainly  from  the  data  given  by  Hall,  et  al.*’~*’.  The  basicity  constants  for  weak 
bases  such  as  dioxane  and  diethyl  ether  were  calculated  by  making  use  of  the 
experimental  formula  given  by  Gordy  and  Stanford  in  their  papers’**’. 


3.  Discussions 

As  in  the  case  of  phenol,  the  fact  that  both  of  the  two  absorption  bands  are 
related  directly  with  the  hydrogen  bonding,  appears  no  doubt  from  their  behaviors 
by  temperature  changes  of  solution  and  concentration  changes  of  proton  acceptor, 
it  further  being  cleared  up  by  considering  the  deformation  of  the  potential  curve 
arisen  from  the  hydrogen  bond  formation*.  It  is  clearly  seen  in  Figs.  1  and  2 
that  both  the  relative  intensities  and  the  wave  number  differences  of  the  two 
absorption  bands  are  closely  connected  with  the  basicity  constants  of  the  acceptor 
molecules.  That  is  to  say,  as  the  base  strength  of  the  acceptor  molecule  increases, 
the  longer  wavelength  component  of  the  two  bands  decreases  in  intensity  with 
the  increase  of  the  shorter  one,  and,  at  least,  in  the  range  of  higher  basic  strengths, 
both  of  these  correlations  reverse.  Such  remarkable  behaviors  of  the  two  bands 
can  be  interpreted  qualitatively,  on  the  basis  of  the  double  minimum  potential 
previously  described,  as  follows. 

Abnormal  changes  of  the  band  intensity  due  to  the  hydrogen  bonding  have 
been  pointed  out  by  several  workers.  Explanation  for  the  line  broadening  might 
be  reduced  to  (1)  the  increment  of  the  effective  charge*’**’  and  (2)  the  increase  in 


3)  N.  F.  Hall:  J.  Am.  Chem.  Soc.  52  (1930)  5115. 

4)  N.  F.  Hall  and  M.  R.  Sprinkle:  ibid.  54  (1932)  3469. 

5)  L.  A.  Flexser,  L.  P.  Hammett  and  A.  Dingwall:  ibid.  57  (1935)  2130.  • 

6)  D.  H.  McDaniel  and  H.  C.  Brown:  ibid.  77  (1956)  3756. 

7)  W.  Gordy  and  S.  C.  Stanford:  J.  Chem.  Phys.  9  (1941)  204. 

8)  W.  Gordy:  ibid.  9  (1941)  215. 

*  The  observed  frequencies  of  the  fundamental  and  its  overtone  bands  of  the  tth  normal 
vibration  in  a  polyatomic  molecule  are,  according  to  Herzberg’s  notation,  expressed  by 

Calculating,  by  using  the  frequencies  previously  measured  by  Badger,  et  al.  in  gaseous  state, 
<u«**  and  xk”  for  the  methanol  hydroxyl  stretching  vibration,  3768  cm~*  and  — 86cm~*  respec¬ 
tively  are  obtained.  The  frequencies  of  the  fundamental,  the  first  and  the  second  overtone 
bands  obtained  in  carbon  tetrachloride,  carbon  disulfide,  dichloroethylene  and  chloroform 
solutions  fit  in  considerably  well  with  the  above  equation,  the  corresponding  potential  constants 
in  carbon  tetrachloride  solution,  for  example,  being  evaluated  as  3726  cm~*  and  —  82cm~‘ 
respectively.  In  the  previous  cited  paper^’,  these  were  assumed  for  phenol. 

9)  H.  Tsubomura:  J.  Chem.  Phys.  23  (1955)  2130;  24  (1956)  927. 

10)  G.  M.  Barrow:  J.  Phys.  Chem.  59  (1955)  1129. 
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E9u  <Pi 


anharmonicity  of  the  potential  curve**' ■**'  arising  from  the  hydrogen  bonding**. 
If  only  the  relative  intensity  of  the  two  absorption  bands  is  concerned,  the  former 
could  be  left  out  of  our  considerations.  The  latter  is  to  be  reduced  eventually  to 
the  modifications  of  the  eigenfunctions  of  ground  and  excited  states  corresponding 
to  the  infrared  absorption  bands  under  consideration,  in  which  mixtures  of  eigen¬ 
functions,  which  belong  to  the  two  potentials,  must  be  superimposed  in  consequence 
of  a  passage  through  the  potential  barrier  or  tunnel  effect. 

Since  asymmetrical  properties  in  double 
minimum  potential  for  weak  base  like  aceto¬ 
phenone  seem  comparatively  large,  relative 
configurations  between  energy  levels  cor¬ 
responding  to  the  two  bands  may  be  shown, 
for  example,  in  Fig.  3,  in  which  the  separa¬ 
tion  of  the  unperturbed  excited  levels  is  not 
very  large  as  described  below.  In  con¬ 
sequence  of  the  perturbation,  the  vibration¬ 
less  eigenfunction  can  be  shown  _  to  be 
and  also  the  perturbed  excited 


Vt 


Fig.  3.  Correlation  of  energy  levels 
for  acetophenone. 


eigenfunctions  to  be  the  following  mixtures  of  the  zeroapproximation  eigenfunctions: 

¥ii=a(Pi+b<Pi 


(1) 


where  the  coefficients  a  and  b  are  determined  as  functions  of  the  perturbation 
energy  W X2  and  the  separation  of  the  unperturbed  levels,  J,  as  follows: 

a  =  (V4\Wt2\^  +  J-‘‘  -  •f/2  v/lJW^pTT*)* 

(2) 

The  perturbed  level  E  is  given  by 

E=}^  iE\+E\)±  A  ( 3 ) 

Then  the  difference  between  the  energies  of  the  resulting  two  states  is  expressed  by 

J£:  =  i/4|IF,*1"-|-J"-  (4) 

From  Eq.  (1),  the  matrix  elements  of  the  transition  moments  corresponding  to 
the  two  absorption  bands  are  calculated  as 

11)  R.  H.  Badger  and  S.  H.  Bauer:  J.  Chem.  Phys.  5  (1937). 

12)  M.  Davis  and  B.  M.  Sutherland:  ibid.  6  (1938)  755. 

13)  H.  L.  Frisch  and  G.  L.  Vidale:  ibid.  25  (1956)  982. 

**  The  spectral  complexity  may  also  be  considered  as  resulting  from  molecules  having 
various  degrees  of  polymerization.  See,  for  example,  Cross,  Burnham  and  Leighton:  J.  Chem. 
Soc.  59  (1937)  1134;  R.  Mecke:  Disc.  Farady  Soc.  9  (1950)  161. 

This  line  broadening  will  be  further  discussed  in  later  paper. 
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Moi  =  jvj*M9'^=a^(Pt*M<f>odv-b^<Pi*MVodr 

(5) 

Moii  =  ^9ii*M9^T=ia^«>i*MVidT+b^^»*Me>^T. 

Under  the  conditions  ft>a,  which  is  obtained  from  Eq.  (2),  and  ^<Pi*M<Paih'> 
^Vx*M(p^r,  which  may  easily  be  understood  from  the  overlapping  between  the 

eigenfunctions,  Mti'^M^n,  which  requires,  as  shown  by  the  experimental  fact,  the 
intensity  of  the  shorter  wavelength  band  to  be  larger  than  that  of  the  longer 
wavelength  band. 

The  basicity  constant  is  considered  as  a  measure  of  proton  attracting  power 
<Broensted  base).  As  this  constant  of  the  proton  acceptor  molecule  increases,  it 
is  to  be  expected  that  both  the  zero-order  frequency  and  the  anharmonicity 
constants  for  the  first  potential  increases,  which  correspond  to  elongation  of  the 
OH  bond  length,  and  that  the  minimum  of  the  second  potential,  on  the  other 
hand,  is  lowered  with  the  increment  in  the  zero-order  frequency.  As  a  consequence, 
the  unperturbed  energy  levels  approach  each  other,  which  would  be  followed  by 
the  reversal  of  their  relative  configurations.  That  is  satisfactorily  inferred  from 
the  fact  that,  with  the  increase  of  the  basicity  constant  of  the  acceptor  molecule, 
both  of  the  two  absorption  bands  shift  toward  the  longer  wavelengths  and  the 
magnitude  of  their  saparation  becomes  minimum  at  an  intermedate  basic  strength, 
as  shown  in  Fig.  2.  The  above-mentioned  facts  are  discussed  somewhat  quantita¬ 
tively  in  the  following. 

Q>nceming  the  band  intensity  and  from  Eq.  (4),  it  is  very  desirable  to  mea¬ 
sure  integrated  extinction  coefficients  rather  than  extinction  coefficients  of  band 
maxima,  the  former  necessitating  a  correction  for  finite  slit  width  to  be  made 
which  is  very  small  in  an  instruments  with  a  large  resolving  power.  If  the  feature 
of  the  observed  absorption  bands  follows  the  Lorentz  curve,  a/fv— ,  in  the  * 
case  in  which  the  difference  in  the  half-widths  seems  unlikely  large,  the  former 
is  possible  to  be  replaced  approximately  by  the  latter,  since  the  integrated  and 
maximum  extinction  coefficients  are  given  by  nalclb  and  ajclb  respectively*^’. 

Neglecting  the  matrix  element  ^<Pt*M*padT  in  Eq.  (5),  the  intensity  ratio  is 
given  by 

=  J)/(V'4  i  W»  '+df-ff)={dE+S)ME-S) .  (6) 

The  separations  of  the  unperturbed  levels,  d,  are  therefore  evaluated  using  the 
values  of  i  and  dE  determined  experimentally  with  the  aid  of  Eq.  (6),  followed  by 

14)  D.  A.  Ramsay:  J.  Am.  Chem.  Soc.  74  (1952)  72. 
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calculation  of  the  corresponding  perturbation  energies  Wu  from  Eq.  (3)  shown  in 
Table  1.  Furthermore,  the  values  of  the  coefficient  b  which  are  obtained  from 
Eq.  (2),  and  the  shifts  of  the  unperturbed  levels  which  are  calculated  by  the  use 

of  Eq.  (5)  from  the  original 
level  corresponding  to  the 
unassociated  band,  Ji®  and 
^2®,  are  also  listed  in  the  same 
table.  In  Fig.  4,  J,  J,®,  Jj® 
and  b  are  plotted  respectively 
against  the  basicity  constants 
of  the  acceptor  molecules. 
As  seen  in  the  figure,  mini¬ 
mum  value  of  J  is  probably 
situated  at  the  basicity  con¬ 
stant  with  which  the  relative 
intensity  of  the  two  absorp¬ 
tion  bands  becomes  equal  to 
one,  which  is  in  good  agree¬ 
ment  with  what  is  likely  to 
be  expected  from  Eqs.  (2)  and 
(3).  In  addition,  for  weak 
base  such  as  nitrobenzene,  the 
basicity  constant  of  which  is 
about  —10,  the  value  of  b  is 
almost  equal  to  one,  which 
results  in  disappearing  of 


100 


Fig.  4.  Correlations  of  J,  Ji®,  J2®  and  b  with  the 
basicity  constants  of  proton  acceptors. 


longer  wavelength  one  of  the  two  absorption  band  agreeing  well  with  the  experi¬ 
mental  fact.  The  perturbation  energies  are  all  approximately  constant,  which 
also  to  be  anticipated  qualitatively,  and  their  values  appear  reasonable  as  compared 
with  the  value  of  157  cm”'  for  the  doublet  splitting  of  2v2  vibration  of  ammonia. 
Moreover,  it  is  of  interest  to  note  that  while  ii®  considerably  varies  with  the 
basicity  constant  of  the  acceptor  molecule,  J2®  remains  approximately  constant. 

Thus  the  remarkable  behaviors  of  the  intensities  and  the  separations  of  the 
two  absorption  bands  related  with  the  changing  of  basic  properties  of  the  proton 
acceptor  molecules  could  be  understood  by  the  deforming  of  the  double  minimum 
potential.  Furthermore,  the  feature  of  the  potential  curve  could  almost  be  deter¬ 
mined  in  combination  with  the  data  which  are  obtained  from  measurements  of 
the  fundam.ental  and  the  first  overtone  regions  and  by  analysis  similar  to  the  above. 
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Finally,  let  us  discuss  the  dependence  of  the  frequency  shift  on  the  hydrogen 
bond  energy.  It  seems  reasonable  to  consider  the  value  of  J  described  above  as  a 
measure  of  the  strength  of  the  hydrogen  bond.  In  fact,  J  monotonically  varies 
with  the  changing  of  the  ba¬ 
sicity  constant  of  the  proton 
acceptor  molecule,  as  seen  in 
Fig.  4.  and,  in  the  case  in 
which  it  can  be  regarded  that 
only  a  single  absorption  band 
results  from  hydrogen  bonding, 
its  shift  from  unassociated  ab- 
soption  band  has  usually  been 
considered  as  a  measure  of  the 
strength  of  the  hydrogen  bond. 

In  Fig.  5,  the  average  value, 

8,  of  the  shifts  of  the  two  ab- 

200  400  600  800 

sorption  bands  from  the  free  j 

OH  band,  which  is  obtained  Fig.  5.  Relation  of  average  values  of  the  shifts 
from  only  the  experimental  the  two  absorption  bands  to  J’s. 

work,  is  plotted  against  J.  The  linear  relation  between  them  as  seen  in  the 
figure  shows  that  the  relative  strength  of  the  hydrogen  bond  can  be  inferred  from 
the  magnitude  of  the  former.  It  is  important  to  note  that,  as  clearly  seen  in 
the  figure,  the  order  of  hydrogen  bonding  strengths  of  the  acceptor  molecules  is 
reverse  to  that  expected  from  the  magnitude  of  their  dipole  moments,  whereas 
it  is  parallel  to  that  expected  from  the  magnitude  of  their  ionization  potentials, 
as  mentioned  in  the  previous  paper". 
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